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We study a three-components fermionic gas loaded in a one-dimensional optical trap at half- 
filling. We find that the system is fully gapped and may order into 8 possible phases: four 2/cf 
atomic density wave and spin-Peierls phases with all possible relative n phases shifts between the 
three species. We find that trionic excitations are unstable toward the decay into pairs of kinks 
carrying a fractional number, Q = 3/2, of atoms. These sesquions eventually condense upon small 
doping and are described by a Luttinger liquid. We finally discuss the phase diagram of a three 
component mixture made of three hyperfine level of 6 Li as a function of magnetic field. 
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A problem analogue to quark confinement in parti- 
cle physics has been recently addressed in systems with 
multi-components fermionic atoms loaded in an optical 
lattice^"^. These studies strongly support the forma- 
tion of "baryonic" states made of bound states of n > 2 
atoms. In one dimension, for example, trionic states, 
made of bound state of three atoms, were predicted to 
be stable at generic densities and sufficiently low temper- 
atures, typically of the order ~ 30 — lOOnKM This re- 
sult opens the exciting possibility to probe in a new con- 
text and in future experiments a simplified version of the 
"quark" confinement phenomenon in quantum chromo- 
dynamics. In all these previous studies the attention has 
been drawn on the formation of baryonic (or molecular) 
states that contains an integer number of atoms. Here 
we shall focus on the intriguing situation where the low- 
energy elementary excitations carry a fractional number 
of atoms. Although it may appears counter-intuitive, 
fractionalization of quantum numbers is a well estab- 
lished phenomenon in condensed matter physics. Cel- 
ebrated examples are fractionally charged excitations in 
the quantum Hall stated and in quasi-one-dimensional 
polymers^. In this work we shall present strong argu- 
ments that fractionalization can also occur in ultra-cold 
atomic physics. We shall give evidences that at densities 
close to half-filling a three-components fermionic mix- 
ture loaded in a one-dimensional optical trap may sup- 
port low-energy excitations carrying a fractional number, 
Q = 3/2, of atoms, the sesquions. 

When loaded in a one-dimensional optical lattice of 
wavelength A, a three-components mixture is well de- 
scribed, away from resonance, by a Hubbard-type hamil- 
tonian of the forrrP^* 

H = ~* £ [4,a c i+l,a + H.C] + Uab W 

i, a i,a<b 

where c\ is the creation operator for a fermionic atom 

of species a = (1,2,3), at site i, and = c\ a Ci, a is 
the local density of the atomic species a. The param- 
eters t and the couplings U a b can be expressed in term 



of the recoil energy, the laser intensity and wavelength 
as well as the s-wave scattering lengths s a b(B) between 
the species. For generic external magnetic fields B, the 
s a b are in general different and so are the couplings U a b- 
Thus the physical symmetry group of ([I]) is U(l) 3 corre- 
sponding to the conservation of the number of atoms of 
each species. Such a small symmetry, which is an essen- 
tial feature of atomic mixtures, make the elucidation of 
the physics associated with ([!]) a difficult task. However, 
as we shall see, much can be said in the weak coupling, 
low-energy, limit. The physics described by ([I]) strongly 
depends on the density of atoms p — pi — p2 = p~3- 
Away from half-filling, i.e. when p ^ 1/2, it was shown 
in Ref. 11 that for generic couplings the dominant fluctu- 
ations consists into massless 2kp Atomic Density Waves 
(ADW) and massless trionic excitations carrying total 
atomic number Q = 3. At half filling, when p = 1/2, the 
physics is, as we shall see, radically different. 

Effective Low Energy Hamiltonian. The low energy ef- 
fective theory associated with the Hubbard Hamiltonian 
([I]) can be derived, as usual, from the linearization at the 
two Fermi points ±kp of the dispersion relation of free 
three-component fermions: 

Ci, a ~ *aR ^ + *aL e~ ikFX a = (1, 2, 3), (2) 

where x = iao, ao = A/2 is the lattice spacing, A the laser 
wavelength, and kp = np/ao is the Fermi wave- vector. 
Finally p = 1/2 is the density per species. In the weak 
coupling limit \U a b\/t << 1 the effective hamiltonian as- 
sociated with (fTl) is found to be: 

H = H +H 1 (3) 

with 

«0 = -™FB*U*aH-44*J ( 4 ) 

a 

and 

= { flab ha R hbL + Kb CbB^bL + Kb J lbR J abL) 

a<b 

+ (R^L) (5) 
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FIG. 1: Generalized atomic density wave and spin-Peierls 
phases. Atomic species are labelled 1,2 and 3. Among the 8 
possible phases j4+ ++ , , A+_ + , ^4^ ++ four are repre- 
sented. In the Air r_ phase the species 1 and 3 are in phase 

and out of phase with the species 2. 

where 

Kr(L) = ^ aR{L) ^aR(L) , llm{L) = *L(L) ^bR(L) , (6) 

and 

^lbR(L) = ^Ir(L)^Ir(L)- (7) 

The fermi velocity is given by up = 2t, \i a b = U a b, 
= ±U a f, and we have omitted a term ^2 a<b (h a RhbR + 
habhbL) that account for a non- uniform velocity renor- 
malization. In the absence of the A + -term we recover 
the hamiltonian studied in Ref. [11] in the p ^ 1/2 
case which involves the 9 U(S)\ R (i l ) currents J\\r(l) = 
( h aR(L),ll bR{L) ,I abR{L) )- The A+-term is a 4/c F con- 
tribution of the density-density interaction among the 
species and is present only at half-filling. It involves the 
6 currents J X r(l) = { J lbR(L)> J abR{L)) that creat e or de- 
stroy pairs of atoms. Together with the f/(3)|fl(n cur- 
rents they generate the 15 SO(6)\ R ^ currents Jr[l) = 
J\\R(L) © J±r(l)- The non interacting part of 
is that of relativistic free fermions and has the maxi- 
mally available symmetry SO(6)r (g> SO(6)l generated 
by the Jr{l) ' s - The interaction hamiltonian breaks the 
later symmetry down to U(l) R ® t^(l)||diag correspond- 
ing to the conservation of the number of atoms of each 
species a. Though the lattice hamiltonian ([T]) depends 
on three couplings U ab the hamiltonian ^ is the most 
general hamiltonian for the three species problem with 
an U(l) 3 symmetry and one has to consider the role of 
the 15 couplings (A^ h , p ab ) that encode all possible com- 
peting orders. Which one is likely to be stabilized in the 
low energy limit depends on the asymptotic behavior of 
the Renormalization Group (RG) flow. 

Renormalization Group Analysis. We have obtained 
the one-loop RG equations associated with They 
will be given elsewhere^ and we shall only present in 
the following our results. Due of the lack of symme- 
try in the problem, and consequently of the relatively 
large number of independent couplings (A^ fc , p a b), it may 



seems an akward task to draw any general conclusions 
on the phase diagram associated with Fortunately, 
it has been recognized^ that due to the importance of 
strong quantum fluctuations in one dimensional systems 
the symmetry at level of the lattice spacing ao is likely to 
be enlarged at low energies which thus considerably sim- 
plify the problem. Phrased in the RG language such a 
Dynamically Symmetry Enlargement (DSE) correspond 
to a situation where the hamiltonian ^ is attracted un- 
der the RG flow toward an effective hamiltonian H* with 
a higher symmetry. As shown in Ref.HTlthe possible DSE 
fixed points H* depend only on the symmetry breaking 
pattern described by ([3| . In the present case we find that 
for generic initial conditions of the RG flow, (A* fc ,^ a fc), 
the low-energy physics associated with Q is described 
by one of the fixed points hamiltonians: 

Ke 2 e 3 = -™J2(<R d **aR-*l L d**a L ) 

a 

± G fce a (¥ aR * aLT *l L * a R?J (8) 

where v is a renormalized velocity, e a = ±1, and G is 
some positive coupling. As Q is invariant under the 
simultaneous change e a — > — e a , there are 8 independent 

fixed points, with (£16263) = (+++), ( — h+), (H h) and 

(++ — ), that describe phases, labelled Af ie2e , with quali- 
tativally different physical properties. The phases Af ie2E3 
and Aj ie2€3 , are generalized 2kp, Spin-Peierls (SP) and 
Atomic Density Wave (ADW) phases. The e a account for 
all possible relative 7r-phase shifts between the species a. 
A pictorial representation of the ground states is pre- 
sented in Figjl] The corresponding lattice order param- 
eters can be readily obtained from the structure of the 
interacting part of ^ and are given by: 

a + l£2£3 = E ^^(4^+1 + hx) (9) 

a 
a 

In each phase the ground state is doubly degenerated 
and when < C^ e2£3 >^ there is spontaneous sym- 
metry breaking of translational invariance by one lattice 
site. As a consequence we expect kinks (or solitonic) ex- 
citations that interpolate between the two ground states 
to be present in the spectrum. As we shall see these have 
fractional quantum numbers. 

Spectrum and Fractionalization. Remarkably enough, 
the different hamiltonians (JsJ) can be brought to the same 
form by mean of duality transformations 1 -: 

Ht ie2e3 (G,*)=UX ++ (G,cjt ie2e3 m, (10) 

where the duality transformations to act only on the right- 
moving fermions as: 

^e 2 e 3 (*afl)=e^ (lTl) e * Qfl . (11) 
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We therefore find that the elucidation of the low- 
energy physics described by the fixed points hamiltonians 
Ht ie2C3 stem from the knowledge of those of H+ ++ . The 
latter hamiltonian has an enlarged 50(6) symmetry gen- 
erated by the J A = J dx(j£ + A = (1, 15), 
and is that of the 50(6) Gross-Neveu (GN) model. 
The other fixed points hamiltonians ([8| has a dual ex- 
tended symmetry 50(6) generated by the dual currents 
J A = / dx{u£ ltata (j£) + j£). Fortunately, the SO (6) 
GN model is integrable 15 so that its spectrum is ex- 
actly known and by duality the one in the other phases 
^4^e 2e3 - In all cases it consists into two set of four kinks, 
5^ =0 i2 3i °f mass mj ~ ie - */* 7 (U being a character- 
istic energy scale), belonging to the two spinorial repre- 
sentations of SO (6) (or 50(6)) and six real (Majorana) 
fermions £p=i . e of mass mp = y/2 ms transforming 
according to the vectorial representation of 50(6) (or 
50(6)). Though their wave functions are different in 
all the phases Af ie2e3 , these particles are described by 
the same quantum numbers since the duality transfor- 
mations (11) do not affect the 3 conserved charges (or 



Cartan generators): 



Qa = J dx (*l R * aR + *l L * aL ), a = (1, 2, 



3). (12) 



These are nothing but the total number of atoms 
of a given species a and the particles of the spec- 
trum are labelled by the set of quantum numbers 
<?2> (fe)' The kinks quantum numbers are fractional 
: 5± = ±1/2(1, 1,1), 5± = ±1/2(1,-1,-1),^ = 
±l/2(-l,l,-l) and 5^ = ±l/2(-l, -1, 1). The 
fermions are bound states of two kinks and have the 
same quantum numbers as the original lattice fermions: 
^ = (±1,0,0),(0,±1,0) and (0,0, ±1). There are no 
other stable particles. In particular there are no stable 
trions in contrast with what happens at incommensu- 
rate fillings. Trionic excitations, = cjc^Cg, have total 
atomic number Q — 3, where 



Q = <7l +<72 +<?3, 



(13) 



and are unstable toward the decay into elementary kinks 
or fermions. A trion has quantum numbers (1, 1, 1) and 
the most energetically favorable process is — > Sq Sq 
so that one may think of the kink Sq as "half a trion. 
As it has total atomic number Q — 3/2 one may call 
it a sesquion. The existence of these fractional kinks as 
the lowest energy excitations in generic three-components 
ultra-cold atomic systems is an unexpected and non- 
trivial finding and constitute one of the main results of 
the present work. It is therefore worth discussing the 
stability of the above excitations. Indeed, as noticed in 
both Refs. (|16I17|) . the DSE is only approximate and we 
expect residual symmetry breaking operators to survive 
even in the low-energy limit. As the 50(6) GN particles 
are labeled by the conserved quantum numbers ( 12 ) as- 
sociated with the t/(l) 3 symmetry of the problem, small 



residual anisotropy will result only into a small splitting 
of the particle spectrum. For large enough anisotropy 
and/or strong couplings it may eventually happens that 
the above description of the spectrum breaks down. We 
expect however that the DSE description of the hamilto- 
nian ([I]) (and hence the stability of the sesquions) holds 
in a relativly large portion of the phase diagram. Indeed 
the accuracy of the DSE description has been checked 
numerically in the different context of the 51/(4) Hub- 
bard model at half-filling where the adiabatic continuity 
of the 50(8) GN spectrum in this case has been explicitly 
observed for small enough interactions^. 

Doping. To model small doping we consider adding a 
chemical potential term Hq — —fiQ (we consider hole 
doping with /i > 0). As Q is invariant under the duality 
transformation ( 11 ) it is sufficient to consider doping the 



50(6) Gross-Neveu model. The chemical potential term 
breaks the 50(6) symmetry but since [H+ ++ ,Q] = 
doping does not spoil integrability and the following pic- 
ture emerges. At non zero /i, the particle spectrum is 
splitted according to the values of Q. A particle with 
mass m and atomic number Q will lower its energy to 
m — nQ. When its energy becomes negative the ground 
state start to fill with these particles. When p > 2ms /3 
the first particle that start filling the ground state is the 
sesquion Sq of mass ms and Q — 3/2. As /i is in- 
creased further other particles would like to enter the 
ground state like other members of the kinks multiplets 
with Q = 1/2 or the fermions with Q = 1. However 
the increase of the chemical potential is counteracted by 
the repu lsion felt by the kinks and the fermions to the 
sesquion d 19 l 20 [ As a result for p > 2ms/ 3 the ground 
state is only filled by sesquions which become massless 
excitations. The effective theory describing these mass- 
less fractional Q = 3/2 excitations is a Luttinger liquid 
with a stiffness K. At these dopings the kinks 5„_j 2 3 
and the fermions £p remain massive, with renormalized 
masses. Both the renormalized masses and the stiffness 
K could be in principle computed from the Bethe ansatz 
solution in a similar way as done in Ref. 1191 At large 
doping, i.e. when p >> ms, the Akp term A + in (pM) 
decouples and one recovers the physic described above 
for the generic filling case with massless Q = 3 trionic 
excitations. We may therefore expect that below some 
critical value of the density p < p c sesquions get con- 
fined into trions. The elucidation of the nature as well as 
the location of such a confinement / deconfinement transi- 
tion goes beyond the scope of this work and will studied 
elsewhere^. 

Three-Species Problem and Experiments. The 
phase diagram in the three-dimensional space 
{Ui2/t,U23/t,U~3i/t) is rich and complex, revealing 
the delicate balance between the different competing 
orders. We find that among the 8 possible phases only 
5 are stabilized: the four ADW phases A~ C2C3 and 
the uniform SP phase We show in Fig. 2. two 

projections of the phase diagram in the (U23/t,U 3 i/t) 
plane for typical value of U^/t = ±0.5. Though it is 
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a) b) 



FIG. 2: 2D projections of the weak coupling phase diagram 
of the anisotropic Hubbard model for fixed values of Unit: 
a) Uiz/t = -0.5 and b) Uvtjt = 0.5. The SP phase A+ ++ 
occurs in the frustrated region (i.e. when (IT a< b U a t) > 0) 
around the isotropic points \U a b\ = U . 

difficult to draw any general quantitative picture of the 
phase diagram we observe that: i) in the unfrustrated 
regions, where (II a< & U a b) < 0, the ADW phases that are 
stabilized are the one that minimize the density-density 
potential in ([I]) and ii) that in the frustrated regions, 
when (n a< f, U a b) > 0, the kinetic energy term play 
an important role when all couplings are of the same 
order of magnitude. In particular in the vicinity of the 
isotropic rays \U a b\ — U a uniform SP phase ^4+++ 
is likely to be stabilized. Eventually the above SP 
phase is destabilized in favor of various ADW phases 
for sufficiently large anisotropies. In experiments once 
the optical lattice parameters and the density are fixed 
the only control parameter is the external magnetic 
field B. The phase diagram as a function of the 
magnetic field B is a line in the three-dimensional space 
(U 12 {B)/t,U 23 (B)/t,U 3 i(B)/t) which dependence on B 
essentially depends on the mixture through the s-wave 
scattering lengths. Taking for example^ a mixture made 
of a balanced population of three hyperfine states of 



6 Li atoms, \F,m F ) = |1) = |l/2, 1/2), |2) - 1 1/2, -1/2), 
and 1 3) = |3/2, — 3/2), with typical optical lattice 
parameters^ and a laser wavelength A = lpm, we find 
a weak coupling regime where a non trivial SP phase 
may be observed. At half-filling using our one loop RG 
equations the following phase diagram as a function of 
the magnetic field emerges. For small fields B < B cl , 
a uniform ADW phase A+ ++ is stabilized while at 
larger fields B > B c2 an ADW phase AZ ++ , where the 
species labeled 1 is in phase opposition with species 2 
and 3, shows off. This match the result of Ref. [TT] for 
p 7^ 1/2 where at these filings ADW phases of the same 
type with quasi-long range order were predicted. The 
essential difference with the above case is that when 
p = 1/2 an intermediate uniform SP phase is 
locked in the region B c2 < B < B c \. Within the one 
loop accuracy we find B c2 ~ 560G and B c \ ~ 540G. 
This is an interesting result from the experimental point 
of view since in both ^4+ ++ and AZ+j. phases we expect 
that effect of the three-body losses^ will be considerably 
reduced. The knowledge of actual values of the binding 
energies of the sesquions and hence of the temperature 
scale below which these phases could be stabilized 
call for an alternative approach such like numerical 
calculations^. To summarize we have shown that in 
the vicinity of half-filling fractional excitations carrying 
Q = 3/2 atoms, the sesquions, are the relevant low 
energy excitations in a generic three-components Fermi 
mixture. These sesquions are likely to get confined into 
Q = 3 trionic excitations when one moves sufficiently 
far away from half-filling. We therefore expect that 
both the confined (trionic) and unconfined (sesquionic) 
phases could be probed in future experiments. 

We thank E. Boulat, S. Capponi, V. Dubois, P. Le- 
cheminant and M. Najac for discussions and encourage- 
ments. 
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